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Abstract 

We consider a one dimensional interacting particle system which describes the 
effective interface dynamics of the two dimensional Toom model at low noise. We 
prove a number of basic properties of this model. First we consider the dynamics 
on a finite interval [1,1V) and bound the mixing time from above by 2 N. Then we 
consider the model defined on the integers. Because the the interaction range of the 
rates and the jump sizes can be arbitrarily large, this is a non-Feller process. As 
a result, we can define the process starting from product Bernoulli measures with 
density p G (0,1), but not from arbitrary measures. We show that the only possible 
invariant measures are those product Bernoulli measures, under a modest technical 
condition. We further show that the unique stationary measure on [0, oo) converges 
to i.i.d. Bernoulli variables when viewed far from 0. 


1 Introduction 


In this paper, we consider an interesting interacting particle system originally introduced 
in [5] to describe the effective dynamics of the interface between two phases in Toorn’s 1 
Model (also known as the North-East, or North-East-Center, model) in the limit of weak 
noise. We recall here (see [13] for more details) that Toom’s model is a discrete time 
probabilistic cellular automaton on Z 2 in which the spin configurations at £ { — 1, l} z2 
are updated according to the rule 




' sign (a t (i,j + 1) + a t (i + 1 ,j) + a t (i,j )) with probability 1-p-q 
< +1 with probability p 
— 1 with probability q. 


The parameters p, q represent noise in the update scheme. It is remarkable, and important 
for what follows, that, for p and q small enough, the system has two stationary states, 
one with mostly +l’s and the other with mostly — l’s. 

1 Pronounce Toom with a long o, not with the English pronunciation of oo 
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One may impose an interface between these two phases by setting the model up in the 
third quadrant of 1? and fixing boundary conditions for a t (j, 0) = +1 and = —1 

for all j < 0 and for all t. 

lip = q = 0, all “up-right” paths from (—oo, —oo) to (0,0) define stable configurations 
(with + above and — below the path) for the deterministic dynamics and one may ask 
how these interfaces fluctuate for p,q 0 but small. Heuristically, one may expect that 
flips off the line separating the + and — regime die out quickly, and the dynamics of the 
line is governed by flips on it. For example, a spin-flip at a vertex immediate to the left 
of long vertical segment of the interface will generally cause further spin-flips at vertices 
adjacent to the segment of interface below. The net effect on the interface is to shift 
the part of the segment below that point 1 unit to the left. Encoding vertical edges of 
the line by +1 and horizontal by —1, the authors of [5] arrive at the following effective 
description of the dynamics at weak noise. First of all, the up-right paths are encoded 
by spin configurations a := (cr(x)) 3;e N G { — 1,1} N (1 corresponds to a vertical segment 
of the interface and —1 to a horizontal segment). Second, the dynamics on up-right 
paths is described by a continuous time Markov chain. Each 1 particle is equipped with 
an exponential rate A_|_ clock, and a —1 particle with a rate A_ clock. We will assume 
throughout that A + , A_ >0 and that A + + A_ = 1, the latter condition simply fixes of 
unit of time. When the clock rings for a particle of fixed sign, the particle exchanges 
positions with the first particle to its right of opposite sign. From now on we will refer 
to these dynamics as the Toom Interface. We will not return to the two dimensional 
dynamics. Here and below, we will denote this process by at := (a t (x)) x ^. 

One remarkable feature of this model is that its restriction to the first N vertices 

{± 1 }“. N1 is itself a Markov chain; the dynamics is the same unless a clock rings for 

a spin in the last block of constant sign in {1,..., N}. For updates of spins in the last 
block, the dynamics reduces to single vertex spin flips. In the language of [5], there are no 

finite-size effects. It is easy to check that the restricted chain is irreducible on {±1}T. 

and hence has a unique stationary measure 7 Tjy. The sequence of measures {ttn)ngn is 
consistent, and this in turn implies that the full chain has a unique invariant measure on 
{±1} N , 7TOO, which restricts to ttn on {±1}{ 1 ’-"’ JV J'. 

Very little is understood rigorously regarding the behavior of either 71^ or the process 
at, though the papers [3, 5, 6] contain a number of interesting conjectures, heuristics and 
numerics. The first paper on the subject, [5], studied the Markov chain defined above as 
a model describing fluctuations via kinetic roughening, the height function h x (at) being 
defined by h x (a t ) = Yli=i a t{i)■ The striking observation there is that if A + = 1/2, the 
statistical properties of the model cannot be in the class governed by the conventional 
KPZ equation: in this case the process h x (at ) is distributionally invariant under global 
spin flip. 

One way to understand this at a heuristic level is to follow the work of Kardar, Parisi 
and Zhang and guess the behavior of h(at) in the appropriate scaling limit. The process 
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should satisfy the SPDE 


d t h = kA h + W(t, x ) + aiyh) 2 + 6(Vh) 3 


where W is a space-time white noise and the last set of terms make explicit the possible 
dependence on the gradient of h. If A + / A_, one concludes that only the quadratic 
term is relevant using scaling theory [7]. However, if A + = h and —h are identically 
distributed, which forces a = 0 in any putative scaling limit. The extent to which 
the third order term is relevant is an intriguing open question. It is marginal in the 
renormalization group sense, and as such [5, 6] argue against its appearance for the 
scaling limits of microscopic models. The situation here is analogous to the expected 
relationship between the scaling limit of the Ising model in 4 dimensions and the putative 
4>\ field theory. 

The simplest manifestation of the above discussion appears in the study of the vari¬ 
ance, under n^, of the sum of the first L spins. Numerics, Renormalization group calcu¬ 
lations and heuristics [5, 12] suggest that 


Var,,- 


oo 



L 2/ 3 if A + + i, 
L 1 / 2 log 1 / 4 L if A+ 


l 
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It might help the reader to consider what this implies about the correlations in this 
model. If Var <C L then the model must exhibit strong negative correlations to cancel 
the contribution to the variance coming from the term ^ a 2 = L. 

With this background in mind, our paper constitutes the first rigorous analysis of the 
Toom interface, though the results fall short of answering the most intriguing questions 
raised in [5], e.g. the above hypothesis on the variance. (The paper [2] analyzes a similar, 
but different model). Let us now present our main findings. We recall that the total 
variation distance between two measures /i, v on a finite sample space tt is defined as 

IIm — u \\ ■= l IM a ) - KcOI 


Abusing notation slightly, we also use at to denote the restriction of the chain to {±1}^ 
and let fi) denote the distribution of at when starting from the initial configuration 
£ € Recall that the mixing time of at is defined as 

Tmix(-V) := inf 

Our first result is as follows. 

Theorem 1.1. For all N € N, 

Tmix(A0 < 2 N 
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A potentially surprising property of the Toom interface model is that it can be de¬ 
fined on the whole of Z. In this case the Bernoulli i.i.d. measures are invariant to the 
dynamics. In other words, the phenomenon of unusually small variance seems to dis¬ 
appear (notwithstanding that one expects, as in ASEP, to recover small variances when 
studying certain dynamic observables such as the current across an appropriately chosen 
space-time characteristic). We wish to understand this disparity better. We will do it in 
two different directions. 

The first direction is to study the behavior of n <*, in the bulk, far to the right of 0. 
Is it Bernoulli? Note that on Z all i.i.d. product Bernoulli measures, with any density 
p, are invariant to the dynamics. On N far from the boundary, the prospective Bernoulli 
measure is fixed by the condition [a x \ = p (this being dependent on A + and A_). 

Formally, let t x be the translation by x i.e. for any spin configuration er, with domain 
D C Z, let r x a denote the spin configuration with domain D + x defined by ( T x a)(y) = 
a(y — x). Denote the induced map on the space of probability measures by r*. Studying 
the behavior of 7roo to the far right is thus studying limfc_ > ._ 00 t^'k 00 . 

Theorem 1.2. Consider (Tf.ir 00 )k£-'N as a sequence of probability measures on {±1} Z . 
Then this sequence converges weakly, as k —>• —oo, to the i.i.d. Bernoidli measure Ber p 
with 

) ’ = £ 

The second direction is to ask: are there any measures on Z invariant to the dynamics 
other than the i.i.d. Bernoulli measures? We show that none exist, under some conditions 
which promise that information does not flow too fast from — oo. While we failed to 
construct an “exotic” (i.e. non-i.i.d.) invariant measure, we have no good reason to 
conjecture such an example does not exist, it seems a condition on the flow really is 
necessary. As the specific conditions we use are somewhat lengthy to state, we defer the 
statement of this result, Theorem 2.6, to the next section. 

1.1 Proof ideas 

The main tool that we employ is a coupling. Let a 1 and a 2 be two starting configurations. 
We wish to construct a coupling of the Toom processes starting from a 1 which makes them 
attempt to become similar with time. We perform this coupling as follows. We start with 
independent Poisson clocks (one for each vertex) each with rate 1. Suppose there is a 
Poisson arrival at time t and at a site x. We examine cr^(x) and cr 2 (x). If u}(x) = of (a;) 
we want the particles at x to move together. To obtain the proper particle clocks we have 
to reduce the rate, so we throw a coin with probability 1 /2 and make them both walk if 
it succeeds (for this informal discussion we assume A + = A_ = l /2, the A + 7^ A_ case is 
similar). If (t}(x) 7 ^ a 2 (x) then we again throw a coin with probability 1 / 2 : if it succeeds 
we make cr 1 walk, and if it fails, we make a 2 walk. It is easy to check that both < 7 ) are 
Toom processes, so this is indeed a coupling. 
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Let us examine discrepancies i.e. x such that cr l {x) 7 ^ <j 2 (x). A Poisson arrival at 
x will force a 1 (a:) = cr 2 {x) after it, but a discrepancy might form somewhere to the right 
of x, call this site y. We say that the discrepancy at x moved to y. Our discrepancy 
might also move because of a Poisson arrival before x, but the key point is that it in all 
cases it moves to the right. The important points regarding discrepancy dynamics are as 
follows: discrepancies are never created. When they move, they only move to the right. 
They may annihilate each other, but only by collisions between opposite types; e.g. a “+ 
discrepancy” (a discrepancy where (j 1 (x) = 1 and er 2 (x) = — 1 ) hits a discrepancy”. 
Let us note here that this coupling is attractive: If Uq < Uq pointwisely, then they remain 
so for all future time. 

Let us sketch how the coupling gives our results. 

Sketch of a proof of Theorem 1.1. Recall that we want to show that the mixing time on 
a finite interval of length N is bounded by 2N . We couple two processes on this finite 
interval with arbitrary starting configurations and examine the discrepancies. They move 
right with speed bigger or equal to 1/2 and fall off the right edge. By time 2 N they are 
all gone, and the configurations are the same. This is well-known to imply a mixing time 
bound. □ 

Sketch of a proof of Theorem 1.2. Recall that we wish to show that the Toorn process 
on N, examined at x, is approximately Bernoulli. We couple the half-line process to the 
full-line process with the correct p. In this cases the coupling may create discrepancies. 
For, a Poisson arrivals at some non-positive x which moves a particle to some y € N does 
not have a half-line process counterpart (note that whether a discrepancy is created or 
not depends also on whether a 1 (y) agrees with u 2 (y) prior to the arrival). We wish to 
analyze the flow of discrepancies across a half-space deep in the bulk. We therefore move 
to a version where the coupling is stationary too (we already have that both coupled 
processes are stationary, but the coupling is not necessarily so). This is done using a 
more-or-less standard limit process. We then examine the rate at which discrepancies 
flow past a point x € N (denote it by j x ). We show that j x is a decreasing function 
of x , with j x — j x+ 1 being exactly the rate of annihilations at x. To prove the result, 
it is enough to show that show j x tends to 0. Since the rate of annihilations must be 
small, the only way that j x cannot tend to 0 is if there are, with positive density, long 
stretches (in space or in time) of discrepancies of a single sign. However discrepancies can 
only annihilate in pairs, so long stretches of discrepancies of the same sign correspond to 
periods of time in which the signed sum of discrepancies across 0 is large. Finally, we 
show that the latter cannot happen often enough to support a non-zero limit for j x . □ 

Sketch of a proof of Theorem 2.6. The theorem will state that the only stationary mea¬ 
sures on Z are Bernoulli. To show this, we start with a stationary measure /j and couple 
it to all Bernoulli processes at once (themselves coupled so that for every p < q the 
Bernoulli-*/ process is pointwise larger than the Bernoulli-/? process). To couple more 
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than two Toom processes at once, just do as follows: once the site to move is selected, 
throw a random coin, it it falls on heads move all Is, and if it falls on tails move all —Is. 
As in the previous proof sketch, we construct a version where the coupling itself is also 
stationary. We then show that there cannot be any annihilations in the coupling, as the 
flow of discrepancies is stationary, and annihilation would cause the set of discrepancies 
to decrease with time (we need here that the flow is finite and space-bounded, which 
induces some conditions on our measure p ). This means that, compared to any of the 
Bernoulli-p measure coupled to it, it is either pointwise bigger than it everywhere, or 
pointwise smaller. There is, thus, a critical p (possibly random) such that p is Bernoulli- 
p (perhaps except at one point). From here it is not difficult to conclude that p is a 
mixture of Bernoulli measures. □ 

Of the three, the most accessible is the proof of Theorem 1.1 appearing in § 2.3. § 3 
is devoted to a proof of Theorem 2.6 while § 4 gives a proof of Theorem 1.2. Except 
for some notation set out at the beginning of § 3 these latter two sections may be read 
independently of one another. § 5 contains a number of lemmas used in both § 3 and § 4, 
notably the existence of a stationary coupling. 

Let us also mention a second paper in preparation [4], In that paper, we prove 
various functional central limit theorems for additive functionals of local observables, 
local currents, tagged particles and the like. Combining the results of that paper with 
the present paper, we are in fact able to derive the bound 

L 

Var Woo ( a *) ~ L ■ 

X=1 

Going beyond this bound probably requires a new idea beyond the technology developed 
here and in [4]. 
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2 The Main Coupling and Dynamics on {±1} Z 

The heuristic given in the run-up to Theorem 1.2 presupposes that the dynamics may 
actually be defined on Z. This is a nontrivial issue as as the process does not have a finite 
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interaction range - arbitrarily distant parts of the configuration on the negative axis 
can influence the local jump rate - and hence the standard Hille-Yosida construction, as 
outlined e.g. in [9], is not applicable. As far as we know, only a few non-Feller interacting 
particle systems have been constructed, most of them relying on a monotonicity property 
that is missing here, see e.g. [11, 8]. But beyond applicability of standard tools, there 
are serious issues of existence and uniqueness. The process is not defined starting from 
arbitrary starting conditions in any reasonable sense - how would one go about defining 
it if the starting conditions are, say, + on the entire negative line? Further, even if for a 
given starting configuration and collection of Poisson arrivals there exists a version of the 
process which is defined for all time, it is not clear that such a version is necessarily unique. 

For example, suppose the configuration is ••• + -(-1- H-- and that there is a 

Poisson arrival at —2 n at time 1/n, for all n (with —2 n the first in the block of two signs). 
Then forcing the spin at —4n to jump 2 units and all those at —4n+2 to jump one unit is a 
legal solution, but so is its opposite. We have no example of an initial configuration where 
such non-uniqueness occurs with positive probability and constructing such an example 
should be interesting. To circumvent this we start by requiring uniqueness, which is 
encoded in Definition 2.1 below. But first some preliminaries. 

Rather than thinking of A-|--Poisson clocks as being attached to particles, we will 
consider a sequence of i.i.d. rate one Poisson point processes (N x (t)) x associated with 
vertices x G Z. Besides these Poisson point processes, the sample space on which our 
coupling is defined supports a two dimensional array of of i.i.d. uniform [0,1] variables 
{U x ,j)x£Z,jeN- Let (17, Pdenote a probability space which supports all these vari¬ 
ables. Define the filtration of sigma algebras {^t)te R+ on by 

&t = cr (N x (s) : s <t; U x>k ■ k < N x (t)). 

Finally, let D = D([0,oo) —> {±1} Z ) be the space of cadlag functions from [0,oo) to 
{±1} Z where {±1} Z is equipped with the usual product topology. We equip D with the 
Skorokhod topology. 

Definition 2.1. A T-process is a pair (p,F) where p is a probability measure on {±1} Z 
and F is a function from {±1} Z x 17 to D which is Borel measurable. F need only be 
defined pxP-almost everywhere. We require F(r],uj)(0) = r) and further that {F(r],u)(t) : 
t G [0, T]} be APt- measurable. 

Alternatively, a T-process is a D-valued random variable a such that for some (p, F) 
as above, P(er G E) = (p x P)(F” 1 (£')). 

We denote F t {ipu) = F(rj,uj)(t ) and similarly a t is the random variable on {±1} Z 
given by a at time t. 

Nothing has yet been formalized regarding Toorn processes in this definition, so we 
just called it a “T-process” in anticipation of the Toorn model, which enters in the next 
definition. We remark that a clearly determines p as erg = p, and determines F p x 
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P-almost everywhere, which is enough, since F is anyway defined only p x P-almost 
everywhere. 

Definition 2.2. For a G D, u G i?, t £ [0, oo) and x < y in Z, we say that a Toom jump 
is required for a at (t, x, y ) if the following occurs: 

1. There is a Poisson arrival at (t,x) i.e. N(x,t) = N(x,t~) + 1. 

2. If a x (t~) = 1 then U x , N r t ) is required to be less than \ + , otherwise it is required to 
be bigger than A + . 

3. a x (t~) = a x+ i(t~) = ■■■ = <x.y_i(f _ ) = —a y (t~). 

A Toom process on T, is a T-process (/a, F ) such that p x P -almost surely, F(rj,uj) has the 
following properties 

1. Ft(r},ui)(x), considered as a function oft, has only finitely many jumps in any finite 
interval, for any x. 

2. If a Toom jump is required for F(ri, u ;) at (t, x, y), then there are spin flips at x and 
y at time t. Otherwise there are no jumps at time t. 

A Toom process on S for some S C Z is the same, except that in Clause 2 we also require 
x€S. 

We are now in a position to define our coupling, which is simply using the same u> to 
run a number of different Toom processes. Formally, 

Definition 2.3. Let {(fl l ,F' 1 ) : i G 1} be two or more Toom processes (not necessarily on 
the same subset ofL). When we discuss a “coupling of the (p l ,F l ) started from /. i” we 
mean the following: p, is assumed to be a measure on Ilie/ {±1} Z whose marginals are 
the ji l . The coupling is then the collection of the D-valued random variables a 1 given by 

a 1 = F l (if,!jj) a 1 : (Y\{±l} z ^) x 12 ^ D. 

iei 

The coupling with independent starting positions is the object given when p is taken 
to be pi. 

The reader is urged to verify that this definition corresponds to the verbal definition 
given in § 1.1 (this has to do with the precise way in which U enters in the definition). 

2.1 Statement of Theorem 2.6 

Theorem 2.6 is about stationary Toom processes, so let us start by defining those. The 
definition uses the natural time shifts on fl. Denote them by St, i.e. St as the map which 
takes N x (■) i —y N x (- + t) — N x (t ) and similarly for the JJ. We can now define 



Definition 2.4. A T-process (n, F) is stationary if 

1. F preserves fi, i.e. for any t € (0, oo) and any E C {±1} A Borel, 

J 1 E(Ft{v,u))dn{y)dP(u) = p,(E). 

2. F forms a semi group i.e. F s (F t (rpu), Stco) = F s+t (ri,uj) for all t,s € (0, oo), /ixP- 
almost everywhere. 

It is called “stationary on S ” for some S C Z if we only require clause 1 to hold for events 
E which depend only on S, formally if rj e E,rj\s = rj'\s => r f £ E. 

We need one more technical condition. 

Definition 2.5. A T-process (n, F ) is called regular if for every t one can write F t as the 
/i x P -limit in measure of functions F t L such that F^(-,u>) is continuous (as a function 
from {± 1 } Z to itself), for almost all uj. 

We have no example of a Toom process which is not regular. Heuristically, construct¬ 
ing a non-regular example seems a similar challenge to constructing an example of a 
non-unique Toom process (recall the discussion on page 7). We will not do it here, but it 
is possible to formulate very mild conditions of “no flow of information from infinity at 
finite time” which would ensure that a process is regular. We are not very happy about 
this condition, but at least, as will be shown in § 2.2, it is very easy to check in concrete 
cases. It will be used only once, in Lemma 3.6 below. 

Finally, let l y (a ) and r y (a) denote the cardinality of the maximal block of spins with 
the same sign to the left of y , starting from y — 1 and respectively to the right of y, 
starting from y + 1 , in particular l y ,r y > 1 . 

Theorem 2.6. Let a be a stationary, regular Toom process on 7L. If the integrability 
condition 


supE[(/ a ,) 1 +e (cr 0 )] < oo, 

holds for some e > 0, then no is a mixture of product Bernoulli measures. 

By making additional assumptions on a, e.g. spatial translation invariance, one may 
assume weaker moment conditions. However to keep the presentation streamlined, we’ll 
stick with the above hypotheses on all measures encountered below. 

2.2 Construction of standard processes 

There is no problem in constructing a Toom process on a finite interval or on a cycle: 
There is only one choice of F which is to go over the Poisson arrivals in the interval 
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one-by-one and apply the rules stated above. Thus F is well-defined unless there are two 
Poisson arrivals at the same time or infinitely many Poisson arrivals in a finite interval 
of time. Since both have probability 0, our F is defined P-almost everywhere, which, as 
already mentioned, is enough. 

Turning to a discussion of the evolution on viewed as a factor of {±1} N , 

we restrict our attention to initial configurations having both infinitely many +’s and 
infinitely many — ’s in N (a property which is preserved for all time), Ft is a Markov chain 

on . N }. It is straightforward to check that this is an irreducible Markov chain on 

a finite state space, therefore it has a unique invariant measure. (To show irreducibility, 
note that to get from a configuration rj to another, ?/, first make the Poisson clocks of all 
— sites in i] ring from right to left, getting to the all + configuration. Then have all — sites 
in rf ring from left to right, getting to rf). With this initial measure, (n,F) is stationary 
on {1,..., N} in the sense of Definition 2.4 (the fact that F preserves ji and satisfies the 
and the semigroup property follows from a direct check of the pertinent definitions). The 
process is regular because it is in fact itself continuous for all oj for which it is defined. 

Extending the definition to a Toom process on N is almost as easy. Let F N be the 
Toom process on {1,..., N} defined in the previous paragraph (recall that the definition 
assumed that the initial configuration has infinitely many + and infinitely many — in N). 
Then F N are consistent in the sense that for x £ {1,..., N}, F t N (rj, u)(x) = F t M (rj, u)(x), 
for all M > N, all t, all r] and almost all c o. Hence their stationary measures — denote 
them by f.i N — are also consistent. Taking a limit gives a measure ji°° and a function 
F°° and it is straightforward to check that (^ 00 ,i' 10 °) is a Toom process on N if only /r°° 
gives zero measure to configurations with a tail of a unique sign. 

Lemma 2.7. jF° gives zero measure to configurations with a tail of a unique sign. 

Proof. Fix N and let M S> N. Examine the event in \.lm that cro(N) = ao(N + 1) = 
• • • = (To (M) and assume for concreteness that the common value is +. The process exits 
this state with rate at least A +(M — N ), as any Poisson arrival in {N,.... ,M} with an 
appropriate U will exit this state. However, it is not difficult to check that returning 
to this state requires at least one Poisson arrival in {1,..., IV — 1} or in M. So its 
rate is bounded by N. We get that the probability of this event is bounded above by 
N/(X + (M — N )). Taking M —> oo shows that the probability of <Jq(N) = a(N + 1) = • • • 
in fj ,oo is zero. As N was arbitrary, the lemma is proved. □ 

The Toom process just described is regular because it is a limit of F t N —>• Ft also 
/i°° x P almost surely. It is also easy to conclude from the stationarity of (/i' v , F N ) that 
(/r 00 ,! 700 ) is stationary. 

The last result we wish to show here is that i.i.d. Bernoulli-p processes have a cor¬ 
responding F which makes them into a stationary, regular Toom process on Z. This 
is a known folk result, but it seems it never appeared in writing so we put it here for 
completeness. 
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Lemma 2.8. There is an F : {±1} & x fl —> D such that for all p £ (0,1) the couple 
(Ber p , F) is a stationary, regular Toom process. 

Let us isolate the first step of the proof as a separate claim. 

Lemma 2.9. To prove Lemma 2.8 it is enough to construct such an F which has the 
required properties only for t < e for some fixed e. 

Proof. Exchange e and 2e, and call the input of the lemma G , i.e. G : {±1} Z x 17 —}D 
and it has the required properties only for t < 2e (t + s < 2e for the property that it 
forms a semigroup). We form F by repeatedly applying G i.e. 

F t (rj,u) = G t (rj,u) t<e 

F t (ri,uj) = F t ^ e (F e (r},u),S e u}) t>e 

where S e is the time shifts on 17, as in the previous section. Note that we are using here 
stationarity: iy_ e is defined only p x P-almost everywhere, so the expression only makes 
sense because the couple (F e (r],uj), S e uj) has /ixPas its law. 

It is easy to check that F preserves p and that it is a Toom process on Z. To check 
that F forms a semigroup, assume first that t < e and get 

F s {F t (r],u),S t u) = F s _ e (F e (F t (r!,u;),S t uj),S e+t u) 

'---' 

= F s _ e (F t+e (p,Uj),St + eUj) 

= F s _ e (F t (F € (rj,u),S e io),S t+e u ) 

= F s+ t- e (F e (r],uj),S € uj) 

= F s+t (r),u). 

The first equality follows from opening the outer F by its inductive definition, the second 
is obtained by applying the semigroup property to the inner term (marked by a brace), 
which is allowed since t + e < 2e. The third is reopening in the opposite order of t and e. 
The fourth is by assuming the semigroup property has been proved inductively for s — e 
and the fifth is again the definition of F. This shows the case t < e by induction on s. 
Concluding the case of general t is similar and we will skip it. 

Finally we need to show that F is regular. We show that by induction on t so we 
will assume it has already been proved for t and will demonstrate it up to t + e (our 
assumption on G is the induction base). Fix 5 > 0 and let F t L satisfy that 

p x P({(» 7 ,w) : d(F t (ri,uj),F t L (ri,u})) < >1-^5. 

Since F t L (- , ui) is defined on a compact space, it has a modulus of continuity (which might 
depend on u). Take such a modulus which holds for all but probability, i.e. a 7 which 
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satisfies 


P({w : d(r],rf) < 7 =► \F t L {p,u) - F t L (i/,w)| < ^5}) > 1 - ^<5. 

Finally pick M such that 

P x P({fa,w) : d(Gf(ri,uj),G e (ri,uj )) < 7}) > 1 - ^<5. 

Thus F^{G^(rj,uj), S e uj) is a 5-approximation of F t (G e (r),oj), S e u), which, by the defini- 
tion of F, is the same as Ft+ e . And of course, it is continuous for almost all u>. This 
finishes the lemma. □ 

Proof of Lemma 2.8. We will construct F by taking the limit of finite systems with pe¬ 
riodic boundary conditions. Let us define the system with periodic boundary conditions 
formally, even though there are no surprises. We define p = p\ L "* to be the following 
process on given as a function of Poisson arrivals: suppose we have an arrival 

at time t and position x. Let y be the cyclically first point right of x having opposite 
sign, i.e. pt-{y mod 2 L) = —pt-(x) and pt-(z mod 2 L) = pt-(x ) for all x < z < y (here 
and below, y mod 2 L is the element of (— L, L\ congruent to y modulo 2 L). Now define 

' Pt-(y) z = x 
Pt(z) = p t -(x) z = y 

^pt-(z) otherwise. 

(for completeness let us stipulate that if po is the configuration with all + or all — then 
pt = po for all t). A simple check shows that for every k 6 {0,...,2L}, the measure 
which is uniform over configurations with exactly k + signs is stationary. Hence so are 
the Bernoulli-p measures. Denote by F L the map {±1} Z x f2 —>• D which realizes this 
process on (-L, L] and freezes the configuration outside the interval. 

We will now show that Fj J (rp ui) converges Ber p x P-almost surely for all p as L —> 00. 
This will construct F, show that it is regular and that it preserves Ber p . By the previous 
lemma, it is enough to show this claim only up to some small fixed time e. We will choose 
e later. 

Compare therefore F l {j],oj ) and F L+1 (p,uj) restricted to some small spatial interval, 
say [— K, K]. The following is a sufficient (if far from necessary) condition for them to 
be equal on [— K, K}: there is an x G (— L,—K ) such that no particle passed over x in 
either F L or in F L+l . If we show such x exist, the lemma will be proved. We call such 
x regeneration points. 

To show the existence of regeneration points, examine the following parameter 

—K—l 

h ■= 

i=l—L 
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(recall that r, is the size of the block of spins to the right of i). Since pt is Bernoulli, 
we see that for some C (depending on A + ,p only), It < CL with probability bigger than 
1 — e~ cL (from here and on we call this “with exponentially large probability”, even if c 
might depend on e). We now claim that also 

pf max I t < CL] > 1 — e~ cL . (1) 

Vte[o,e] ~ > 

To see this, switch to discrete time (i.e. examine a process on a cycle that simply moves 
one particle at every step) and notice that, again with exponentially large probability, 
the continuous time interval [0, e] corresponds to no more than L steps of the discrete 
time process. At each step there is exponentially large probability to have It < CL and 
summing over these L times the probability remains exponentially large. This shows (1). 
Denote the event of (1) by 

Order the Poisson arrivals in the interval (— L, —K ) according their time of arrival. 
Assume the k th Poisson arrival caused the particles at x and y to exchange. Define 
Ak = \x — y | and let A be the sum of Ak for all arrivals in [0, e]. Under 'Ll we have 

mi) < c 

and hence A is, with exponentially large probability, bounded by a sum of 2 eL independent 
variables with bounded 4 th moment. This shows that 

P(A > C\eL) < 

But A bounds the number of non-regeneration points (for F L ): Since the same estimate 
holds also for L + 1, we may fix e = l/4Ci and get that with probability larger than 
1 — CL~ 2 a regeneration point may be found. By Borel-Cantelli this means that there 
exists an Lq such that for all L > Lq a regeneration point exists. As discussed above, 
this proves that F is a regular Toorn process that preserves Ber p . 

To finish the lemma we need to show the semigroup property for F. up to time e. Fix 
some K sufficiently large such that F L has a regeneration point in [. K , 0) with probability 
larger than 1 — 5 independently of L (as long as L > I\). Since F L —> F we see that F has 
a regeneration point in \K, 0) with probability larger than 1 — 5. Since 5 was arbitrary, 
we get that F has infinitely many regeneration points. But this shows the claim because 
after a regeneration point we can calculate F as if it were a finite Toorn process. Since 
this satisfies the semigroup property, so does F, and the lemma is proved. □ 

2.3 Two simple applications 

In this section we give two applications of the coupling. The first, Theorem 1.1, was 
already stated in the introduction. The second provides exponential decay of correlations 
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in time. Here is the precise statement: 

Theorem 2.10. Let f,g be local functions with Ber p (f) = Ber p (<7) = 0 and Ber p (/ 4 ) = 
Ber p ((/ 4 ) = 1. Then 

E B er p (/M<7(^)) < Ce r ~ ct , (2) 

with r the length of the smallest interval containing both Supp / and Suppg, and C,c 
only dependent on A±. 

From now on, constants C. c throughout the paper will be allowed to depend on p, X± 
without further mention. 

Proof of Theorem 1.1. Fix <f £ {d=l}^ 1, ’ ,iV ^ to be arbitrary, and let if £ {dzl }^ 1 ’-’ N } 
be distributed according to the stationary measure /cy. Start two Toorn processes on 
{1,..., N} from (f and if and couple them as above (call the resulting processes of and 
of respectively). If we show that at some time T that P(ct^ = of) > \ then T upper- 
bounds the mixing time by definition. 

As explained in the introduction, a crucial property of our coupling is that it pushes 
discrepancies to the right. Let us formalize this statement. Let t\ be the first arrival 
on site 1 and for all j > 2 let tj be the first arrival on site j after Tj-\. Because we 
are looking at the process with a wall to the left of site 1, once t\ occurs, the value of 
of( 1) = o^(l) for all t > T\. By induction, the same is true for all {(j. Tj) : j < N}. 
The theorem is proved by observing that ry is the time it takes for the IV’th arrival of a 
Poisson point process which has rate 1. □ 

Let us remark that there is also a coupling-less version of the argument. Indeed, 
examining only one process, after n the value of <7*(1) is independent of the initial con¬ 
figuration, and similarly for all Tj. Thus ry is a forget time and this is equivalent to the 
mixing time, see [10]. 

Proof of Theorem 2.10. Take two local functions /, g and, for concreteness, say that 
Supp / U Supp g C [0, r]. We define the measure u on o = (a 1 , a 2 ) £ ({±1} A ) as 
follows: 

1. Both o 1 and o 2 are Ber p -distributed. 

2. For x < 0, o l (x) = o 2 (x). 

3. {o l (x) : x > 0} is independent of o 2 and idem with 1 ■£>■ 2. 

Let ot = (of of) be the coupling of Toorn processes started from v (recall Definition 2.3). 

Let X(ot) denote the position of the “left-most discrepancy” of the configuration o tl 
i.e. X(ot) '■= min{x : <r 4 (x) f o 2 (x)}. We are interested in X(ot) for the following 
reason. Since the support of / lies in [0, oo), f(of) is independent of Oq and therefore 
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of of Vt G R + . Moreover, if X(at) > r then <7(of) = g(cr't ). Since f(a q) and <7 (of) are 
independent, we have 




M/^oM )1{^K) > ?•}) +E i/ (/(o-(i)£f(cr t 1 )l{X(o- t ) 


E v(/(^o)g(|> r}) +E i/ (/(a < i) 5 ((cr t 1 )l{X(o- t ) 

e ^(/(^o)5-K 2 )) - e ^(/K)5K 2 )1{^(^) < r}) 

+ E u (f(a^)g{al)l{X{a t ) < r}). 


< r}) 

< r}) 


The hrst term is in the final equality is simply 0 (by the independence explained above), 
so we get 


\^v{f{vl)g(.°}))\ < 2 E Ber p (/ K) 4 ) 7 E Ber p (g K) ^ ' F(X(a t ) < r ) 1 / 2 

< 2P (X(a t ) < r) 1/2 

where we used Cauchy-Schwarz twice for each term, the invariance of Ber p and finally 
that Ber p (/ 4 ) = Ber p (g 4 ) = 1. 

Finally, observe that X(a s ) is naturally coupled to a Poisson process N(s) which has 
rate min(A±) so that 

X(a t ) > N(t) 

The theorem now follows from a large deviation estimate on the Poisson process N(t). □ 


2.4 Generators, Local Rates and Derivatives 

One unfortunate consequence of the fact that the Toom interface is non-Feller is that 
intuitive reasoning involving generators and locally defined rates need to be checked 
rigorously. Here we state two results which justify our pervasive use of this language 
throughout the text, as we will almost always assume Condition A below. Because we 
feel these statements are more of technical rather than actual interest, their proofs are 
relegated to the end of the paper (see § 5.3) 

We will define our operators on the space of continuous functions from {±1} 5 , S C Z, 
to R. We start with the flip operator F x at site x, defined by 

F./W-/M ■>*<»> = v * x 

[-<%) y = x. 

For a finite subset S C Z and a G {—1, l} 5 , we have the indicators 

Xs = xW(x) = d(x) Vx G S] 
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and whenever a is all 1 or all —1, then we simply write Xs an( i Xs- We a ^ so nee d the 
associated projections - 

Psf{?) = Xs( cr )/( £J )- 

Then the generator of the process is formally defined as 

* = E- 2 ’™ = E^+^-nA + A - p [ ;,„-n p v + )( F * F » -1). 

x<y x<y 

We call this definition formal because Jz?/ is in general not continuous, due to the infinite 
sum over x. We need a condition on moments that will almost always be assumed: 
Condition A. We say that a Toorn process a satisfies the uniform Condition A if for 
some r/ > 0 

supE(^ +T? (fJt)) < oo. 

t,X 

We say that it satisfies the local condition A if 

supE (ll +r, (a t )) < oo Mx. 

t 

If you are reading the online version and ever forget what is Condition A, clicking the 
letter A should send you to the definition. 

Lemma 2.11. Let a be a Toom process satisfying the local Condition A. Let f be a local 
function. Then for almost all ao we have that t E(/(er t ) | ao) is differentiable in t, the 
sum defining (Jzf/)(ao) converges, and 

^E(/(cr t ) | £7- 0 ) = (-S? f)(a 0 ) (3) 

Further, we have an averaged version, 

^~E(/(cp)) = E((Jz?/)(cto)) . (4) 

The same result also holds for a coupling of finitely many Toom processes a 1 : if all of 
the a 1 satisfy the local Condition A, then the differentiability of local functions also holds 
for the joint process, with a corresponding formal generator. The proof of this claim 
follows the proof of Lemma 2.11 verbatim. 

Throughout the paper we will also use statements that are slight generalizations of 
the above. The most prominent example comes in Definition 3.2 and Lemma 3.3 where 
we handle “instantaneous rates” that should be justified similarly to the rates of change 
of f{(Jt) above. We will henceforth use reasoning involving rates and currents without 
explicit justification. 
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3 Invariant Measures on Z 

In this section we investigate invariant measures for the Toorn interface on Z, using the 
coupling from § 2. Throughout the section we assume that a = (cr 1 , cr 2 ) is a pair of Toorn 
processes on Z coupled as in Definition 2.3. 

We formally define discrepancies. Let 

D' ? = := {x £ Z : cr^x) = r], cr 2 (x) = —r/} r/ = ±1, 

D = D ff :=D+UD; = {ieZ: a\x) + cr 2 (x)} 

Rather than focusing on the discrepancies themselves, it is useful to restrict attention 
to the study of gaps between consecutive discrepancies of type (+, —) and type (—,+), 
i.e. elements of D + and D~, respectively. Let us, arbitrarily, refer to the first type 
discrepancies as having signature/sign + and the second type of discrepancies as having 
signature —. To keep track of “interfaces” between the two types of discrepancy, let, for 
x € D^, 


6(x) = inf{y > x : y £ D v } 


should such a y exist and set b{x) = oo otherwise. The set of interface discrepancies is 
then 

B = B(a) := {iED: b(x) < oo and (x, 6(x)) D D = 0}. 

The following lemma is a main ingredient of Theorem 2.6, and the place where the 
integrability condition is used. 

Lemma 3.1. Assume (a 1 , a 2 ) is a stationary coupling of Toom processes which satisfies 


supE( m&x.{l x (<j 1 ),l x (<j 2 )}m.m{r x (a 1 )r x (a 2 )}') 1+e < oo. 

X 


Then 

P(cr 1 < cr 2 or a 2 < a 1 ) = 1 

Here and below, “<” stands for pointwise inequality for all t £ [0, oo) and all x £ Z. 
We stress that the phrase “(cr 1 ,;? 2 ) is a stationary coupling” means not only that each 
one is a stationary Toom process, but also that the coupling is stationary. Formally, let 
v be the measure on ({±1} Z ) and F l : {±1} Z —>■ D be the maps that define a i.e 

p(<7 £ E) = J t E (F 1 (r] 1 ,u),F 2 (r] 2 ,uj)) dn(r] 1 ,r] 2 ) dF(u). 

Stationarity here refers to the requirement that (F^rj 1 , uj), F 2 {j] 2 , cj)) is distributed ac¬ 
cording to v for all t. 
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Proof of Lemma 3.1. Denote for brevity 

maxlr = max{l x (a 1 ),l x (a 2 )} minr x = min {r x (a 1 ),r x (a 2 )}. 

In general, our goal is to show that E[|2?ij = 0 for v the measure on ({±1} 2 ) 2 induced 
by do. Let us first sketch the proof in case u is translation invariant. We fix an interval 
I. By stationarity, we are tempted to write 


0 = dt& u [\B n 7|] < E[maxl m ; n /] — min(A±)E[|.B 1 n J|], (5) 


where B 1 C B is the set of interface discrepancies that can be annihilated in one step. 
The inequality in (5) is due to the observations that the first term bounds from above the 
flow rate of of discrepancies from (—oo,min/ — 1] into / while the second term bounds 
from below the annihilation rate inside /. By hypothesis, sup x E[maxl x ] < oo so that 
these two inequalities together imply that Efl-B 1 n /|] is uniformly bounded in I. Under 
the assumption of translation invariance of v, this implies E^B 1 !] = 0. This argument 
can then be iterated (considering the discrepancies that can be promoted into B 1 in one 
step, etc. ). Eventually one then concludes that E[|7>|] = 0. Note that this argument 
is slightly formal due to the fact that|7> fl I\ is not a local function and therefore the 
inequality in (5) would need additional justification. Instead of remedying this directly, 
we proceed to Lemma 3.4, which we need for the general proof of Lemma 3.1 and which 
immediately settles the translation invariant case. We will need the following definition 
together with a justifying lemma, here and in §4. 

Definition 3.2. Let a 1 and o 2 be two coupled Toom processes, and let x £ Z. We define 
j x (a o) to be the infinitesimal rate at which a discrepancy jumps from (—oo,x) to [x, oo) 
(we include also in j x also the case that the discrepancy annihilates in [x,oo) in the same 
step). Formally, for y < x, i £ {1,2} and t > 0 we let J y .t,i be the event that a 1 had a 
Toom jump at time t at position y, and that 1) for all z £ [y,x) <r*_(z) = cr*_(y) and 2) 
for some z £ [y,x), a }-(z) / a 2 (z). Then we define 

j x {a o) = lim -P(3y < x,t < e and i £ {1,2} such that J y ,t,i)- 


We define a x (ao) as the infinitesimal rate of annihilation of discrepancies atx. Formally, 

1 , ■ . 
a x (ao) = lim -P(3 y < x, t < e and i £ {1,2} s.t. J v ti and al_(y ) = of_ l (x) a l .(x)). 
e—>0 e r T 1 

Finally define 

j x = E (jx(cro)) a x = E(a x (cr 0 )). 

These rates are well-defined, as we state now. 
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Lemma 3.3. If a 1 and a 2 are two Toom processes satisfying the local Condition A, 
then for any coupling of a 1 and a 2 , j x (a) and a x (a) are well-defined almost surely and 
their averages a x ,j x are finite. Further, j x and a x (which have been defined above as the 
expected quenched rates) are equal to the corresponding annealed rates, i.e. 

j x = lim-E (if {discrepancy jumps from (—oo,x) to [x,oo) before time e}) 

£—>■0 e 

a x = lim -E iff {annihilations at x before time e}). 
e—>0 e 

Proof. These claims are variations on Lemma 2.11 as a x (cr),j x (cr) are defined as condi¬ 
tional time-derivatives. The only difference with Lemma 2.11 is that these derivatives 
are not naturally equal to E cr («Sf / /) with 2z?' the formal generator of the coupled process, 
but the proof of Lemma 2.11 (given below) applies here as well. For concreteness, let us 
give the expression for a x (a) 

2 

a y( a ) = ^ 

x<y t]=± i= 1 

where V x is the event that g = a 2 (x) = a l (x + 1) = • • • = a l (y — 1) = —a l (y) while 
a 3- * satisfies that a 3 ~ l (y) = g and for some z £ [x,y) we have a 3-l ( z ) = —g. □ 

We now pick up the threads of the proof of Lemma 3.1, starting with 

Lemma 3.4. Let v be the initial measure of a stationary coupling and assume the local 
Condition A. Then 

lim sup E [\B n /|] = 0. (6) 

Proof. Given a discrete interval [x, y], define the event 

E x , y {a) := {x £ B(a), y = b(x)} 

In words, E x>y denotes the event that there is a boundary discrepancy at x and that the 
first discrepancy to its right occurs at y. Recall that a y is the rate of annihilations at y 
and that j x is the rate of discrepancy flow through x defined above. We now claim that 

jx ~ jx+l = a x- (7) 

Indeed, if we denote by H x (t ) the number of discrepancies that crossed x in the time- 
interval [0, t], and A x (t ) the number of annihilations at x in [0, t], then we have 

H x+ i(t) - H x (t) = A x (t) + l(x £ D(cr 0 )) - l(.x £ D(cr t )) 
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Taking the expectation, we can drop the two rightmost terms by stationarity. Dividing 
then by t and taking t — > 0, we obtain (7) by rising Lemma 3.4, more precisely the 
identification of a x as “annealed rates”. In particular, since the j x are finite and the a x 
nonnegative, we have ^ a x < C. 

Then Lemma 3.4 is a consequence of the following claim. 

Lemma 3.5. With v as in Lemma 3-4, and for any x and y, 

y 

22 a z > c(y — x)v(E X) y). (8) 

Z=X 

Let us finish the proof of Lemma 3.4 and then attend to Lemma 3.5. Let / = [y, z} 
be our interval, and let k be some parameter. We get from (8) 

z z+k j 2 

22 u ( x e ■ B ’ b ( x ) - x < k) < 22 f7m ax - C ( 9 ) 

x=y x= y ' 

On the other hand, for any k, 

22 ^ -Bj b(x) — x > k) < 

x=y 


Lemma 3.4 thus follows. □ 

Proof of Lemma 3.5. The main observations we make are that Toom jumps due to ar¬ 
rivals at sites z > y do not harm us — they cannot move the discrepancy at y — while 
arrivals at sites z < x only help us — they can push the discrepancy at x closer to or 
on top of y (such that an annihilation occurs), or can annihilate the discrepancy where 
it stands, but cannot push it beyond y. To arrive at an annihilation event, it suffices to 
have at least y — x Poisson arrivals at the location of the discrepancy which is at x at 
time 0 before any occur between the location of the discrepancy and y. 

To use this we examine the behavior in the time interval [0, 2]. We get 

v(E X) y) = E ^ J l Ex y (<r t ) d t 

< P(3 1 € [0,1] s.t. 1 E x ,y((Tt)) 

< C(y — x )P( ff 1 {annihilations in [x,y\ before time 2}) 

y 

= 2C(y — x) 22 a z 

Z=X 

where the last equality follows from the fact that a z is also the annealed rate of annihila- 
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tions (see the “further” clause of Lemma 3.3). This shows Lemma 3.5 and consequently 
also Lemma 3.4. □ 

We go back to the proof of Lemma 3.1. Let 

B x = {y G B : b(y) < x}, 

and note that 1 {y G B x } is a local function (its support is [y. x\) in contrast to 1 {y G B}, 
so we can use Lemma 2.11. In what follows we fix some x G Z and some h,k G N, and 
we omit all three from the notation to avoid clutter. Set 

1 for y G [x — h, x\, 

0(y) = 0 for y G (—oo, — x — h — k], 

, 1 ~ i for y = -x - h - j, 0 < j < k. 

By stationarity 

X 

0 = StE[ ]T 9(y)t{y G B x } 

y=x—h—k 

Let us examine the events which change the value of the function between [ • ] in time. 
The value is decreased by annihilations and by discrepancies leaving the set B x (an 
annihilation might make the discrepancy just before it become a boundary discrepancy, 
i.e. an element of B x , but since 9 is increasing on (—oo,x] the sum still decreases). The 
latter happens, for a discrepancy at y, when b{y) leaves (oo,x]. We ignore annihilations 
and define 

X := min(A±)l{3y G B x n [x — h + 1, a;] : b(y) can leave (—oo, x] in one step}. 

Next we examine events which cause ^ 9{y)\{y G B x } to increase. Examine one bound¬ 
ary discrepancy y. If r y (o l ) / r y (a 2 ) then, because y is a boundary discrepancy, this 
means y is separated from b(y) by a stretch of equal signs. Thus, a Poisson arrival at y 
either annihilates with the one of opposite type at b(y) or moves one step to the right. 
Similarly, a Poisson arrival in the interval of length [y — maxi y ,y) might cause the dis¬ 
crepancy at y to move 1 step to the right or to be annihilated. The sum in these cases 
increases by no more than 9(y + 1) — 6(y). In the case that r y (a 1 ) = r y (a 2 ) the discrep¬ 
ancy at y might move by this common value. Therefore, in all cases we may bound the 
increase in the sum by 6(y + rnini'y) — 9(y). Defining 

Z(y) := [9(y + minr^) - 6(y)\ maxl y 1 {y G B x }, 
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a convenient way to express the above bounds is to say 

0 < Y E v [Z(y)\-E v X. (10) 

y<x—h 

(note that we use here the existence of the annealed generator, recall (4)). To exploit 
(10), we bound ^ y E [Z(y)\ by splitting the sum over y. For x — h — k < y < x — h, we 
use 1 9(y + minr y ) — 6{y) \ < minr y /k and Holder’s inequality (for any 1 < p < oo with 

l + \ = !) to get 

Y E [\ Z (y)\] ^ (t Y ^[(maxl^minr^n ’ (- Y u( d) G B *) 

y=x—h—k ' y=x—h—k / ^ y=x—h—k 

If we choose q sufficiently close to 1, then the first factor is bounded by C (uniformly 
in h, k, x) by the integrability assumption we placed on v. The second factor decays 
as k —> oo (uniformly in h, x) by Lemma 3.4 and the fact that B x D / C B n I. For 
y < x — h — k, we write 

Y \ maxi x-h-k nihii x - h - k 

y<x—h—k 

since there can be no more than one boundary discrepancy for which 9(y + mmr y ) — d(y) > 
0. We conclude that 

lim sup Y^ nz{y)} = 0 . 

k ^°° y<x 

Combining with (10), we conclude that E„ [X] = 0 (recall that X is nonnegative and 
independent of k). Since this holds for all x,h, it follows that E[|L>|] =0. □ 

At this point we have amassed information about stationary couplings of measures 
but haven’t actually shown any example of such a coupling. The latter issue is the subject 
of the next lemma. 

Lemma 3.6. Let (//,-F 1 ) and, (y 2 ,F 2 ) be two stationary, regular Toom processes. Let v 
be any measure on ({±1} Z ) 2 which has y 1 and y 2 as its marginals. Let u t be the result 
of applying the coupling to v for time t, i.e. 

”t(E) = J t E ((F 1 (jj 1 ,uj)(t),F 2 (r] 2 ,u)(t))dn(r] 1 ,if)dP(uj). 

Then any subsequential weak* -limit of Jp Jq is a stationary coupling of ( y 1 ,F 1 ) and 
(h 2 ,F 2 ). 

As this lemma is technical in nature we postpone its proof to § 5.2. The lemma will 
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also be used in § 4 for constructing a stationary coupling of a process on Z and a process 
on N and the proof in § 5.2 applies to both. Let us mention for now that this is the only 
place in the proof where the Toom process is required to be regular. 

One final lemma is needed before we start with the proof of Theorem 2.6, this time 
unrelated to any coupling. 

Lemma 3.7. Let a be a Toom process on Z (not necessary stationary) satisfying Condi¬ 
tion A. Then with probability 1, if at some time to, the limit 

N h T»2jvTI l{x€[ “ A, ’ ivl: ‘ T ‘« (x, = 1}l 

exists (“(Jo has a density”), then for any t > to we also have that a t has a density, and 
the densities are equal. 

Proof. Fix T and to < T and examine the quantity 

D n := ||{x <E [~N,N] : a T (x ) = 1}| - |{x € [~N,N] : a to (x ) = 1}||. 

Then Dn is the difference between the flows of particles with sign 1 into and out of 
[-N, N] , and since the rates of these flows are bounded by the lengths of the left stretches 
at —N and N we have 


JE-D/v < E / /—Ar(u's) + In(v s ) ds T C 
J to 

(where C depends on the Toom process and on T, but not on N and where the fact that 
the time derivative is bounded by l is again by the annealed generator, see (4)). 

This shows that 2 n+i ^ >n converges in L 1 to 0 for all to < T. To pass to almost 
everywhere convergence we note that if we restrict, say, to squares N = M 2 , then the 
Markov inequality and the Borel-Cantelli lemma together give us that 


1 

2 M 2 + 1 


D M 2 — > 0 


almost surely Vf < T. 


But if 2 N +1 *= [— N,N] : <j t (x ) = 1}| converges as N —> oo on the squares, then by 

monotonicity it also converges with no restriction on the N. Finally, we note that since 
T was arbitrary the claim just proved also holds for all integer T simultaneously, proving 
that 2 N+i ^ n 0 f° r t € [0, oo). □ 

We have gathered all necessary ingredients, and we may now start the 


Proof of Theorem 2.6. Let us fix a stationary measure p, for the Toom interface as in the 
statement of Theorem 2.6. The idea is to construct a coupling P of (u 1 , (cr(j?)) pe [o,i]) such 
that: 
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(a) For p G [0,1], the distribution of a(p) is Ber p . 

(b) The distribution of a 1 is p. 

(c) For any p: P(cr 1 < a(p) or a 1 > (j(p)) = 1. 

(d) If p' > p, then a(p') > a(p) a.s. 

If we then define the random variable 

0 := sup{p G Q : u l > cr(p)} = sup{p : a 1 > a(p)}, (11) 

it is tempting to believe that a 1 = <r(@), and that the distribution of <r(@) is a mixture 
of product Bernoulli’s with mixing measure given by the distribution of 0. To turn this 
into a rigorous proof is a bit delicate. What follows is one possible approach. 

Step 1: Let & C (0,1) be a finite set. Let be a measure on {±1} Z x 

(n ;;e ^{il} Z ) whose first coordinate is distributed like p (this is property (b)) and is 
independent of the others, while its other \&\ coordinates are distributed like {Ber p : 
p G <^} (property (a)) and are coupled to satisfy property (d). Also, let ip denote the 
distribution on {±1} Z x (I~I ; , Gi ^{±l} Z ) obtained by evolving the coupling to time t when 
started from the initial distribution vq. 

We now examine a weak*-limit point of the collection of time averaged 

measures (1/T) / Q T u t d t. The stationarity of Ber p ensures satisfies (a), the stationarity 
of p ensures (b), and of course (d) is also preserved (say because monotonicity is equivalent 
to having only discrepancies of the same type, and this is preserved by the process and 
by taking a weak limit). By Lemma 3.6, is a stationary coupling, so we may apply 
Lemma 3.1 and get property (c). We note that <j(p) is Bernoulli so K[r x (a(p )) k ] < oo for 
all k and x, and ditto for l x . By Holder’s inequality we may thus get from E [^((T 1 ) 1_l_e ] < oo 
(the assumption of the theorem) that E[(maxl J/ minr^) 1 " 1 " 6 ] < oo (the requirement of 
Lemma 3.1). 

To pass from a finite set to an infinite set we only need to make sure that, when 
& C then is a limit of a subsequence of the sequence that was used to define 
This ensures consistency and allows to apply the Kolmogorov extension theorem. We 
may thus get a measure coupling cr 1 to a(p) for a dense set of p in (0,1), say to all of the 
rationals, satisfying (a)-(d). Abusing notations, below when we write pgQwe implicitly 
mean p € Q (~1 (0,1). 

Step 2 : Define the random variable 

0 := sup{p G Q : a 1 > cr(p)j = inf{p : a 1 < cr(p)j 
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The equality follows from monotonicity, and also gives that cr 1 has a density i.e. the limit 

J im wJ—r\i x G i~ N , N } : = 1}| 

AT—>-00 ZiV + 1 

exists, and is equal to 0 since with probability one for all p £ Q the density of a(p) is p. 
Note that the existence of a density is a property of cr 1 irrespective of any coupling of it 
with anything else. 

Step 3: We now wish to claim that the density of a 1 is independent of {<r(p) : p £ Q}. 
For this we need to go through step 1 again. Step 1 starts with the measure v® under 
which cr 1 is independent of {cr(p) : p £ ^}. In particular the density of a 1 is independent 
of the latter collection of random variables. By Lemma 3.7 the density is preserved during 
the time evolution, so 0 exists a.s. under the measures (1/T) v t d t and is independent 
of {a{p) : p £ &*}. 

We need to be careful passing to the limit in T as, in general, weak*-limits do not 
preserve the existence of a density. In this case, however, it is only the coupling that 
changes while the marginal distributions of each coordinate stays fixed. Because 0 can 
be defined purely in terms of cr 1 , for any e > 0 and we can find some local variable E, 
independent of T, which approximates the density of a 1 up to e i.e. P(|.E — 0| > e) < e. 
We get that E is e-independent of {cr(p) : p £ 7^} and this property is preserved in the 
limit. Thus under the density is 2e-independent of {&(p) : p £ As e was arbitrary 
we get the required property for &. This implies directly the property with & replaced 
by Q. 

Step 4: We now show that a 1 = sup{cr(p) : p < 0} under u^. To see this, observe 
that since 0 is independent of {a(p) : p £ [0,1]}, we can fix it in advance. On the other 
hand, for any monotone coupling of Bernoulli processes and any fixed x , the random 
variable X := sup{p : a(p)(x) = 0} is uniformly distributed. Hence for 0 fixed, the 
probability that sup{cr(p) : p < 0} / inf{cr(p) : p > 0} is zero since none of the 
individual coordinates of the monotone coupling can flip at 0. Further, for any value 
of 0 fixed in advance, a := sup{cr(p) : p < 0} is a Bernoulli-© process Since a 1 is 
sandwiched between sup{cr(p) : p < 0} and inf{cr(p) : p > 0}, the foregoing discussion 
implies get that cr 1 = cr(0) and 0 is independent of {cr(p)}. In other words, a 1 is 
distributed as a mixture of Bernoulli processes. This proves the theorem. □ 

Remark 3.8. Surprisingly, perhaps, the standard coupling of Bernoulli processes (i.e. the 
coupling where each site is coupled independently of the others) is not stationary to the 
Toom process (For example, take a finite system with periodic boundary condition. Then 
it is straightforward to write a formula for the number of incoming arrows for any given 
configuration and see that it is not constant). In other words, the process a(p) which we 
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analyzed during the proof is a collection of Bernoulli processes, coupled by a non-standard 
monotone coupling. 


4 Proof of Theorem 1.2 


Recall that we wish to study the process on N away from 0. For this we consider the 
coupling process defined on the configuration space {±1} N x {±1} Z . We will write a = 
(cr 1 ^), cr 2 (x)) for the configuration a € {±1} N x {±1} Z . We use the notation D, D + 
and D~ for discrepancies set out at the beginning of § 3 with the understanding that all 
vertices x < 0 host discrepancies at all times. 

Throughout this section, we let p be the unique solution in (0,1) to the equation 


1 ~P \ 2 = A+ 
P ) A_ 


( 12 ) 


Let v be a probability measure on {±1} N x {±1} Z stationary for the coupling process 
and with respective marginals 7roo,Berp. We know such measures exist by Lemma 3.6, 
see proof in § 5.2. (There is in fact a unique such measure, though we do not need to use 
this explicitly). From the next proposition, Theorem 1.2 follows easily. 


Proposition 4.1. With v as above, 


lim v(x € D) = 0. 

£—>•00 

Before discussing Proposition 4.1 further, let use it to attend to Theorem 1.2. 

Proof of Theorem 1.2. Consider a local function / and recall the shifts t x . By the defi¬ 
nition of D; 


| E »roo[/ or *]-EBer p [/or s ]| < ll/lloo ^ v(y € D). 

yez+Supp/ 

By Proposition 4.1 the RHS tends to 0 as x tends to infinity. But since Ber p is invariant 
under spatial shifts, this implies the push forward of tToo by t x converges weakly to Ber p 
as x tends to infinity. □ 

Proof of Proposition f.l. As in the proof of Lemma 3.4, the crucial quantity that we 
will analyze is the discrepancy flow through a point x, which we will again denote by 
j x . Recall that j x is the averaged infinitesimal rate at which a discrepancy in (— oo,x] 
jumps to (x,oo). Let us note a few simple properties of j. First, u{x G D) < j x because 
whenever there is a discrepancy at x it jumps to (x, oo) with rate at least the rate of the 
Poisson arrivals at x. Next, note that j x is always finite: both cr 1 and a 2 satisfy the local 
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Condition A, a 1 because < x and in particular is finite, while a 2 is a Bernoulli 

process. By Lemma 3.3, this implies that j x is finite. 

Next, again as in the proof of Lemma 3.4, let a x be the rate of annihilations at x, i.e. 
the averaged infinitesimal rate at which a + discrepancy at some y < x moves on a — 
discrepancy at x causing both to disappear, and vice versa. Recall (7) i.e. j x -\ — j x = a x . 
Let k be some fixed parameter and let x € N. One of the following three events must 
occur: 

1. The interval I x := [max(0, x—k 2 ),x\ contains two discrepancies with opposing signs. 

2. The last k discrepancies before x are of the same sign (and case 1 did not happen). 

3. The last k discrepancies before x are not of the same sign (and case 1 did not 
happen) 

We divide the set of discrepancies D and the infinitesimal rates of discrepancy flow j x 
into three parts accordingly, D = D 1 U D 2 U D 2 and j x = j x + j 2 + j' x - Thus, for example, 

D 1 = {x € D : I x n D + / 0, I x f~l D / 0} 

and similarly for the other DL Let us also give an example with j: jf. is the averaged 
infinitesimal rate that discrepancies jump from (—oo,x] to (x,oo) while there are two 
discrepancies of opposite signs in I x (we are not bothered about whether the discrepancy 
that performed the jump is in fact one of those). 

Let us isolate the estimates used to control D,j in two lemmas. We note first of all 
that 

v(x G D®) < j* * = 1,2,3. 

following the same argument we used above to show that v{x € D) < j x . This means 
that our control on D* in these lemmas follows from our control on j l . 

Lemma 4.2. For every k there exists a C*, such that E|D 2 | < Ck- Further, 

OO 

X=1 

The proof (see below) proceeds by observing that a configuration contributing to j x 
has a positive probability to lead to an annihilation. This part of the argument would 
work equally well with any value of p. It is in the next lemma that the specific choice 
(12) becomes crucial. 

Lemma 4.3. There exist numbers e*, | 0 such that ^(|D 2 n[l, ?r]|) < e^n. Further, j x < ek 
for all x > k. 
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We defer the proof of both lemmas and instead first show how they imply Proposi¬ 
tion 4.1. We start with a preliminary claim 

Lemma 4.4. lim n _ ) . 00 ^(|D n [l,n]|) = 0 

Proof. Fix k G N and decompose D into the 3 parts using this k. By the “D parts” of 
Lemmas 4.2 and 4.3 we see that 

^([(D 1 U D 2 ) n [1,n]|) < C k + ne k . 

As for D 3 , the conditions of case 3 imply that the interval [max(0, x — k 2 ),x\ contains 
less than k discrepancies. It follows from this that |D 3 n [l,n]| <n/k + k 2 . This implies 

limsup(l/n)^(|D n [1, n]|) < e k + 1/k. 

n 

Taking k —> oo, the lemma is proved. □ 

To move from the averaged result of Lemma 4.4 to the non-averaged result we argue 
as follows. Let e > 0. First fix some k sufficiently large such that the e k from Lemma 4.3 
satisfies e k < e. By Lemma 4.2 we can choose xq such that for all x > xq one has that 
j], < e. Next, let M be some parameter (to be fixed shortly and depending only on k and 
on A + /A~). By Lemma 4.4 we can further find some x > xq such that 

X 

Y v(v G D) < e/k (13) 

y=x—M 

(we assume here that x > M and, while we are at it, also x > k 2 ). 

We will now show that j x < Ce. Since j], and j 2 are already given to us as satisfying 
such a bound we need only estimate j 3 . For this purpose we note that a 2 is a Bernoulli 
process, so there are constants c,C > 0, such that one cannot create an interval of 
consecutive equal signs larger than Ck before x by changing the sample a 2 at less than 
k deterministic locations, except with probability > 1 — e~ ck . Fix the parameter M from 
the last paragraph to be this Ck and denote this event by G k . Split j 3 = j 3 ' 1 + j 3 ' 2 , where 
j 3 ' 1 denotes the rate at which discrepancies pass from (—oo,x] to (x, oo) by arrivals in 
[.x — M, x\ and j 3 ’ 2 is the remainder. It follows that 

X 

jx' 1 — (M + 1) Y Kye D). (14) 

y=x—M 

because contributions to j 3,1 require at least one discrepancy in [. x—M , x] and also require 
a Poisson arrival in this interval. The factor M + 1 in the estimate comes as an upper 
bound on the latter rate. Using (13) and the fact that M = Ck we get j 3,1 < Ce. 
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To bound j x 2 notice that a pair of particles of opposing signs in both configurations 
acts as a barrier for discrepancies. Hence if an arrival occured at some y and pushed 
a discrepancy beyond x, then either o\ or 02 must have constant sign between x and 
y. Since we are in case 3 there must be less than k discrepancies between x and y (two 
discrepancies of opposing signs are also a barrier), so we have that 02 must have less than 
k particles of some kind. This allows to bound j x 2 by examining <72 only, and 02 is a 
Bernoulli process. The probability that 02 restricted to [y, x\ has less than k particles 
of some kind can be bounded above roughly by Cexp (—ck — c{x — y)/k ) — recall the 
definition of M — and we get 

x—M—1 

Jx 2 < c J2 e~ ck ~ c{x - y)/k < Ce~ ck . 

y=-oo 


Taking everything together we get 


jx <Ce + Ce~ ck . 

and since we assumed k is sufficiently large (depending on e) we may incorporate the 
Ce~ ck into the other term and conclude, as promised, that j x < Ce. 

Proposition 4.1 is now proved because j x are decreasing, by (7). We get that for all 
y > x j y < Ce and since e was arbitrary j y —>■ 0. As already remarked u(x € D) < j x so 
we also get u{x € D) — > 0. □ 

Proof of Lemma 4-%■ An important corollary of (7) is that j x is decreasing and that 
Y^, a x < 00 . The proof of the lemma then follows by comparing j], and a x . j x x is the rate 
at which discrepancies flow through x while there are two discrepancies with opposing 
signs in the interval I x = [max(0, x — k 2 ), x\. This couple of discrepancies acts as a barrier, 
so any arrival that pushed a discrepancy beyond x must be in I x . Hence 

j x < \I X \ ' ^(there are two discrepancies of opposing signs in I x ). 

In the language of Lemma 3.5, this event implies that for some y,z G I x , the event E y z 
occurred. Applying Lemma 3.5 gives 

Jx < k2 u( ' E v^) - fTjE Qw - C ( fc ) aw 

y,z€l x ^ ^ y,zGl x w=y wel x 


as required. 


□ 
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4.1 Proof of Lemma 4.3 

As already mentioned, it is in this part of the proof that the relation between p and A is 
used. We start by indicating the reason for this relation. First denote by H±(t) the total 
number of ± discrepancies which jumped from (—oo,.x] to (x, oo) between time 0 and 
time t. Let H x (t ) = //+(£) + H~(t ) be the total flow and let K x (t ) = Hff(t) — H~(t ) be 
the signed flow through x. We will only need the signed flow at 0 (where it is a function 
of the Bernoulli process alone) and so abbreviate K{t) = Ko(t). 

Lemma 4.5. If ((1 — p)/p) 2 = A + /A_ then E[/f(i)] = 0. 

Proof. This is a straightforward calculation: the infinitesimal rate at which + particle 
enter N is A_|Ji(cr 2 )l{a- 2 (0) = 1} (recall that l\ is the length of the block of identical spins 
from 0 to its left). Since a 2 is a Bernoulli-p process, E[/i(cr 2 ) • l{(cr 2 (0) = 1}] = p 2 . The 
same holds for the — particles and we get E (K{t)) = p 2 \+ — (1 — p) 2 A_. □ 

We next need two preliminary results which are interesting in their own right. It is 
possible to show that K satisfies a functional CLT under proper rescaling, see [4], but for 
our current purposes, the following diffusive bound suffices. Its proof is supplied in § 5.1. 

Lemma 4.6. There exists C > 0 such that for all t € R, 

E Beip [(K(t) - K(0)) 2 ] < Ct. 

We also need a bound on the total flux of discrepancies across a fixed vertex for the 
process on {±1} Z . Set ^V([a,b\) = Ho(b) — Hq(o). 

Lemma 4.7. For sufficiently small 7 , 

EBer p [exp( 7 ^L {I)/\I\)\ < C, for any I. 

In particular 

Peer P (^(I) >N)< C(m)(\I\/N) m (15) 

The proof is also deferred to § 5.1. 

To explain the relation between j 2 and I\ we need the following definition. For a 
partition 7 r of [0, t\ into intervals, let 

Jd(7r) = ^2\AK(I)\l{\AK(I)\>k} 

lEn 

where AK(I) := It (sup/) — /i(inf/) and let 

JH* ( t ) = JF* ( k,t ) = sup JF ( 77 ). 

7T 

With these definitions we can now claim 
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Lemma 4.8. For all x > 1 and t, H x (t.) < x + 

Here Ftf (t) is the number of discrepancies of “type 2” which passed through x until 
time t, using the same classification of discrepancies into 3 types we used on page 27 to 
define j x . 

Proof. The x term in the lemma is a crude bound for “original” discrepancies i.e. for 
discrepancies which existed at time 0 in [l,x], We ignore these discrepancies and order 
the others by their time of crossing x. 

In this proof it will be convenient to think about discrepancies as having a fixed order. 
Examine a discrepancy. When it moves (to the right only!), it can collide with the first 
discrepancy to its right. If the second discrepancy has the opposite sign, the two annihilate 
leaving behind two vertices each having spin agreement. If the two discrepancies have the 
same sign we shall use the interpretation that the first one takes the place of the second 
one, and the second starts moving. This can create a chain reaction, but we note that as 
soon as there is an annihilation, the process stops. In any case, the order of discrepancies 
never changes. 

Let us now use this interpretation to examine a stretch of i type 2 discrepancies 
with the same sign contributing to j x , say with sign + for concreteness. The clock ring 
which precipitates the jump contributing to j x does not have to occur in I x and could 
come from a particle of either sign. However, the chain reaction which occurs can only 
involve the discrepancies with sign +. This follows from our definition; just before the 
jump occurred, the last k discrepancies to the left of x were of the same sign. Since 
discrepancies of different types do not propagate motion when they collide, we conclude 
that the discrepancies of I x must be + discrepancies. 

Hence a stretch of i discrepancies of type 2 which crossed x actually corresponds to 
i + k — 1 discrepancies of sign +, the l discrepancies which crossed as well as another (at 
least) A: — 1 which remain in I x at the time of the last crossing. Let a be the time when 
the first one entered the system and b the time when the last one did. 

We claim that K (6) — K (a) = l+k — 1. This relies on the fact that signed discrepancy 
sum is preserved by our dynamics: discrepancies survive until they are annihilated in pairs 
of + and —. So if we see £ + k — 1 consecutive + discrepancies at some space-time point 
(x,t), we must have started with a signed sum of l + k — 1. (the term “consecutive” might 
be slightly misleading here, since some of them are consecutive in the time of crossing of 
x and others are consecutive in space at a given time, but for the claim on K(b) — K(a) 
this does not matter). 

Finally, different stretches of discrepancies of type 2 crossing x must correspond to 
different time intervals [a, b]. Indeed, examine the last + discrepancy in one stretch and 
the first — discrepancy in the next. Being of different sign they cannot cross one another 
without annihilating, and since we know the + arrived at I x before the —, the — must 
have started after it. 
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Thus the stretches of discrepancies of type 2 crossing x form disjoint collection of 
subintervals of [0. t]. each of which has \AK(I)\ > k. Completing this collection to a 
partition 7 r, the lemma is proved. □ 

Lemma 4.3 will thus be proved once we show 

Lemma 4.9. 

limsup limsup-EBer p [^*(^) ^)] = 0. 

fc—HX> t—»00 t 

Note that we have written EBer p instead of E since K depends only on cr 2 . 

Proof. The proof of this lemma relies on a separation of scales. For each partition 7r of 
[0, t], let us split the sum J^(ir) according to interval sizes: 

JT(7t)= \AK(I)\l{\AK(I)\>k} + J2 |AiL(7)|l{|AiL(/)| > k} (16) 

I(z7r:\I\<k 1 /' 2 

s -v-' s -v-' 

7r) J^ 2 (7 t) 

We will use separate mechanisms to bound each of JP 1 and JP 2 uniformly in 7 r. 

Let us Hrst attend to J^ 1 (7r). The idea here is simple: since a contributing interval 
I is small relative to k , it necessitates too many Poisson arrivals, at least k, in I. This 
is a rare event, and gets exponentially rarer as |/| —> 0. This fact allows us to handle all 
partitions simultaneously via a properly chosen infinite covering of [0, t]. 

For each j £ Z let 


LtPJ 

x j = ^2 Jf{[i2?, (■ i + l)2i))t{jP([i2i, (* + 1)2-7]) > k / 2 }. 

i =0 

Then for any 7 r, exploiting that PP(I) is the total variation of K(I), 


sup J^ 1 ( 77 ) < 2 ^2 x j- 

i<iog 2 (fc 1 / 2 )+i 


Taking expected values on both sides and using stationarity 


£ E Ber p 


sup^ 1 (7r) 


< ]T 2l_JE Ber p K([0, 2' 7 ])l{ t /F ([0, 2-7]) > k}] . 

j<log 2 (fc 1 / 2 ) + l 


(17) 


Using (15), we find that for k sufficiently large, the RHS of (17) is summable and that, 
moreover, it tends to 0 as A: tends to 00 . 

To bound JP 2 {j\) let us introduce a reference partition p = {[ifc 1 / 4 , (i + l)^ 1 / 4 ]} of 
[0,f] (shorten the last interval if necessary). Let p' C p be the collection of intervals which 
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contain an endpoint of some interval from it. Then we can write 

JT 2 (t r) < Ei Air wi+ 2 E jY{I)1{jY{I) > k 1/3 } +2 J2 < k 1/3 } . 



Term III is bounded by AtkT^^kfl 3 since the total number of intervals I £ it with 
|/| > A: 1 / 2 is at most 2tk~ 1 ^ 2 . For the other two terms, we take the expectation and use 
the fact that they no longer depend on 7 r. For II, the bound (15), yields 

E Berp [II] < 2 tk~ 1 ' 4 x C(m)(k 1 ^ 4 /k 1 ^ 3 ) 171 

which vanishes as k —> 00 by choosing m. sufficiently large. For the first term, we argue 

EBerJI] < ^- 1/4 (E B er p [|A'(fc 1/4 ) ~ ^(0) | 2 ]) 4 / 2 < Ctk^ 8 (18) 

where the first inequality is by stationarity and Cauchy-Schwarz, and the second follows 
from the diffusive moment estimate of Lemma 4.6. Hence we have obtained 

limsup(l/f) sup J?T 2 (7r) = 0, 

k—t OO 7T 

Combining with the analogous bound on J^T 1 , the assertion of the lemma follows. □ 

5 Auxiliary results 

5.1 Proofs of lemmas 4.6 and 4.7. 

Proof of Lemma f.6. Let / be the infinitesimal drift of K i.e. f(t) = a(0)\ CT ^' ) li(a) (in 
other words, if cr(0) = 1 then / = A+Zi and otherwise / = —A_/i). Using the Cauchy- 
Schwarz inequality, 

E [K(tf] <2e[(^L-^ /(a s )ds) 2 ] + 2E[(jf /(^)ds) 2 . 

(recall that K and hence / are functions of the Bernoulli process, so all expectations are 
likewise with respect to it). The first term is a martingale, so its expected square is its 
quadratic variation. By direct calculation this is /q A (T (o)Zi(ct), which is bounded by Ct. 
For the second term, we use stationarity to bound it by 

/»oo 

2 1 / |E[/(<t s )/(<7 0 )]| ds. 

Jo 

We shall show E Ber [/(<t s )/(<to)] decays exponentially in s to complete the proof. 
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We want to apply Theorem 2.10. One slight complication is that the function / has 
unbounded support. To handle this issue, let f n be the approximation to / above by 
replacing l\ with rnin{/]. n}. which is a local function. Then, by inspection, 

E[(/ - fn?} < Ce~ cn 

so that, for any n, Cauchy-Schwarz yields 

\nna s )f(a 0 )]-E[fn(a s )fn(ao)}\ < Ce ~ cn / 2 . 

The autocorrelation of f n is handled by Theorem 2.10, with r = n, so that 

\E[f(a t )f(a 0 )]\<C(e- cn / 2 + e n - ct ) 

and, choosing n = ^ct, we get 

|E[/(^)/(a 0 )]|<Ce- ct . 

Lemma 4.6 follows. □ 

Proof of Lemma f.l. Recall that we wish to estimate H = H$(t), the number of par¬ 
ticles entering N in [0, t]. We define H L (t ) to be the number of particles which enter 
from (—L, 0), and note that H L /* H so by monotone convergence E(exp( 7 iL i )) —>• 
E(exp (7 H)). 

Next denote the drift of H L by v = v L (t) i.e. 

:= / min(/i(fj s ), L) ds. 

Jo 

We first Hnd an estimate on E(exp(au)). We Taylor expand the exponential 

OO fc 

^“ v ) = T,'fr E ^ k] ’ 

k =0 A " 

integrate each term and use Holder’s inequality and stationarity to get 

E(v k )< f ... [ t E(h(a sl )---l 1 (a Sk ))ds 1 ---ds k <t k E(l 1 (a 0 ) k ) 

Jo Jo 

so 

E(e Qi; ) < E B er p [e“ Wl ]. (19) 

which is hnite for a < c(p)/t. 
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Next, one verifies by direct computation (omitted here) that, for any (3 > 0, 

Zp(t) ■■= exp {/3H L (t) - (eP - 1 )^} 

is a martingale. To exploit this, fix (3 G R and set a = (1/2) (e 2 ^ — 1) and write 

E Bei> [exp(/3iL i )] = E Bei> [exp (f3H L - au)exp(m;)] . 

Applying the Cauchy-Schwarz inequality and noting that 2 a = e 2/3 — 1 

E Berp [exp (f3H L {t))} < E Berp [Z 2 / 3 (t)] 1 / 2 E Berp [exp( 2 ar;)] 1 ' /2 . 

Since the first term on the RHS evaluates to 1 (being the expectation of a martingale), 
the lemma follows from (19) and E(exp(aLf L )) —> E(exp(aLf)) upon choosing /3 = 7 jt 
and noting that if 7 is sufficiently small then 2 at is small enough to ensure that the right 
hand side of (19) is finite. □ 

5.2 Existence of a stationary coupling 

Our goal in this section is to prove Lemma 3.6. Recall that it stated that if (/it^-F 1 ) 
and (/ 1 2 , F 2 ) are two stationary, regular Toom processes and if v is any coupling of them, 
then any subsequential weak*-limit of y Jq is a stationary coupling of (/r 1 ,F 1 ) and 
(/i 2 , F 2 ). Here v t is the result of applying the coupling to v for time t. As this is the only 
place in the paper where regularity of Toom processes is used, let us recall its definition 
(Definition 2.5 on page 9): a Toom process (^,F) is called regular if Ft can be written 
as the n x P-limit in measure of functions F^ such that F^(-,uj) is continuous for almost 
all uj. 

Proof. As the lemma is used both to construct a coupling between two processes on Z 
(in § 3) and between a process on Z and a process on N (in § 4), let us denote by 5 the 
corresponding state space i.e. ({±1} Z )“ or {±1} N x {±1} Z , as the case may be. Let T n 
be the sequence over which we assume said measures converge weakly* and denote 

1 f Tn 

= — / 17 d t A 00 = lmiA n . 

1 n Jo 

Fix some t. By the definition of regularity, we may write F t ] as a limit in measure of 
functions F f ] ' L such that E t 1,L (-,a;) is continuous for almost all uj . Repeat for F 2 . Let us 
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define the corresponding operators on probability measures on 

P t X{E) = j l B (F i (r/ 1 ,r/ 2 ,w))dA(? 7 1 ,?/ 2 )dP(a;) 

p t L X(E) = j t E (F t L (f? 1 , f? 2 , w)) dA(? 7 1 , r/ 2 ) dP(w) 

where here and below we denote 

F t {ri\ri 2 ,uj) = (F^ (rj 1 , u), F 2 (rj 2 , cv)) 

and ditto for F L . For conciseness, fix t and remove it from the notations Pt and P/\ We 
now fix some / : S —> R and write 

PXoo(f) - Aoo(/) = 

(PXooif) - P L \oo(f)) + (P L Xoo(f) - P L Xn(f )) + (. P L K(f ) - PXn(f)) 

+ (PX n (f) - A.„(/)) + (A n (f) - Aoo(/)) = I + II + • • • + V. 

We will now bound the different terms. We need the following lemma. 

Lemma 5.1. For every continuous f : H —> M and every e > 0 there exists an L such 
that 

j \f(F t L (r ] 1 ,r] 2 ,uj)) - f(F t (r] 1 ,ri 2 ,u))\d\(ri 1 ,ri 2 )d'P(u) < e (20) 

for any A with marginals gf and g 2 and for every t > 0. 

Proof. Let 6 > 0 be some parameter to be fixed later. Let B 1 C H x I? be the set 

B 1 = {(v\v 2 ,u) : | Ft' L (rj 1 ,u)-F?(ri 1 ,u)\ > 5} 

and similarly B 2 with F 2 instead of F 1 . Since F 1,L F 1 in measure, and since B 1 does 
not depend on r/ 2 or on the coupling, we get that for L sufficiently large v x P(F 1 ) < 6, 
and similarly for B 2 . Examine now the integral in (20) and write 


E i.e. 

\/E C E Borel 
\/E C E Borel. 


+ 


+ 


/ B 1 


IB 2 




The first and second terms are each bounded by <5||/||oo since the measures of B l are small. 
The last term is bounded by the modulus of continuity of / i.e. by max{|/(r/) — /(r/)| : 
d(r),r)') < 2d}. If we pick 5 sufficiently small such that the sum of the three terms is 
smaller than e, and the lemma is proved. □ 


We return to bounding the terms I,..., V. To bound terms I and III we use Lemma 5.1 
to get that I, III < e whenever L is sufficiently large (depending on / and e, but inde- 
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pendent of n). Term V converges to 0 as n —> oo because \ n —>• Aqo weakly* and / is 
continuous. For term II we write 


II 


f(Ft'(‘n 1 ,n 2 ,w))d\oo(ri 1 ,Tf 2 ) - / f(F t L (r] 1 ,ri 2 ,uj))dX n {r] 1 ,r] 2 ) J dP(w). 


The functions integrated are continuous (for almost every u), so the inner term converges 
to 0 for almost every ui as n —> oo. It is also bounded since / is bounded. By the bounded 
convergence theorem we get that II —> 0 as n —> oo. 

Finally, term IV is bounded by the observation that PtV s = IJ t+.s■ Let us postpone the 
proof of this fact (which is just playing with definitions) and write 


1 

T'n 





(/) < — max |/| 
-L n 


v s +t ds - 


- > 0 . 

n—> oo 



(/) 


Let us wrap up the calculation. We fix L sufficiently large so that I + III < e uniformly 
in n. We take n —> oo and get |F > A 00 (/) — Aoo(/)| < e. Since e was arbitrary, these are 
actually equal. Since / was an arbitrary continuous function, PtXoo = Aqo. Since t was 
arbitrary, this is the required stationarity. 

We still need to show P t u s = vt+s■ We write 


PtVsU ) = j /(I/(? 7 1 ,? 7 2 ,w))d^(? 7 1 ,r/ 2 )dP(w) 

= f f(Ft(F s (v 1 ,V 2 ,u'),u))dis(r] 1 ,ri 2 )dI > (uj')dP(uj) 

= J f {F t (F s (i] 1 ,u), S s uj)) dz/(r/ 1 , r/ 2 ) dP(cu) 

= J f(F t+s (ri 1 ,r) 2 ,u;))dv(r) 1 ,ri 2 )dI > (u;) = v t+s (f). 


The first equality is the definition of Pf. The second is the definition of v s . The third 
uses that F s is ^-measurable, so we can replace the two independent Poisson processes 
ui and ui' with a single Poisson process, and use co for the first and S s lo for the second 
(recall that S s are the natural time shifts of the Poisson process). The fourth equality is 
the semigroup property for F. □ 


5.3 Proof of Lemma 2.11 

We begin this section with a technical lemma. Let Q(y,t ) be the number of times cr s (y) 
changed in the time interval [0, t] and N y the Poisson process at y. 
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Lemma 5.2. Let a be a Toom process, y G Z and t > 0. T/ien 

EQ(y, t) < f El y (a s ) ds + EiV J/ (i) 

Jo 

(possibly in the sense of oo <oo). 

Proof. Let x < y and let Q(x , y , f) be the number of Poisson arrivals in x in the time 
interval [0, t] that caused the value of a(y) to change (recall from the definition of a Toom 
process that every change in y must correspond to some x < y and some Poisson arrival 
at x at time t such that cr t ~ ( x ) = cr t - (x + 1) = • • • = a t - (y — 1) = — cr t - (y)). Then 

EQ(y, t) = ^2 E Q{x, y, t ) + EN y (t) 

x<y 


and 


El y (a s ) ds = '^ j J E1{<t s (x) = • • • = a s (y - 1)} ds 


where in both cases the change of order of integral, summation and expectation is justified 
by positivity of the integrands. Hence the lemma follows from the next claim. □ 

Claim 5.3. E Q(x, y,t) < E f* 1{cj s (.t) = • • • = a s (y — 1)} ds. 


Proof. Denote the integrand on the right hand side by x( s )- By the definition of a Toom 
process, a s \\ xy \ changes only finitely many times in the time interval [0, t], almost surely. 
Hence x i s Riemann integrable and 


X(s) ds 


[t/e 


= lim e > \{ei). 

e —^0 ^ 


i =0 


Next let t\ < ■ ■ ■ < tk be the Poisson arrivals at x in the time interval [0, t] and let 


B(e) = |i : 3s € (t* - e,U) such that <r a -\[x, y ] + a ^[x, y ]} 


The condition CF s \\ X y\ changes only finitely many times in the time interval [0, t] also 
implies 


lim B(e) = 0 almost surely. 

e —>0 

Hence by dominated convergence (B is bounded by the number of Poisson arrivals at x), 
lim e _>.olE(R(e)) = 0. Defining 

Q = \{i : x(ez) = 1 and 3 tj € [ei,e(i + 1))}| 
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we get Q(x,y,t ) < Q + B (recall that Q requires also the condition a s (y) = —a s (x), 
which is neglected in Q). Since a Toom process has the property that a t - is independent 
of Poisson arrivals at [t, oo) we get 


EQ = E 


"L*/«J 


'L t/e\ 

^2 G [ei,e(i + 1))} 

= E 


i=0 


i=0 


(e + 0(e 2 )). 


Since E(e Yfl'o x{ ei ) ~Jq x( s ) ds ) —> 0 (again dominated convergence), the claim follows. 


□ 


Proof of Lemma 2.11, averaged version (4). Let x < y with y in the support of /, and 
let e > 0 and assume for convenience also e < ^. We wish to define the effect of Toom 
jumps from x to y on /(<j e ). Therefore let t\,... ,tk be the times of the Poisson arrivals 
at x in the time interval [0, e] and define, for x < y, 

k 

D(x,y,e ) = ^2(f(cr ti ) - f(a t -))l{a s (x) = ■■■ = a s (y - 1) ± a B {y)}. 

i =1 

Since, for all z, a Toom process changes at z only finitely many times in any time interval 
and / is local, we get f(cr e ) — /(<7o) = Yl x <y D{x, y, e)- From this we conclude 

E(f(a £ )-f(a 0 )) = ^2nD(x,y,e)) (21) 

x<y 

where the exchange of sum and expectation is justified as follows: Let Q(y, e) be the 
number of times a y changes sign in the time interval [0, e]. Then 

^L>(x,y,e) < 2||/||oo ^ Q(z,e) 

x <y ^eSupp(/) 


which is integrable by Lemma 5.2 and the fact that the sum over the z is finite. Using 
dominated convergence gives (21). 

Moving to the behavior as e —> 0, let us first show that 


d_ 

dt 


E( D(x,y,t )) 


= E(£’x,yf)- Vx,7/ 


t =0 


(where the expectation on the right is with respect to <ro). To see this let E(x,y,e) be 
the event that for some z G [x,y\ and some t £ [0, e], <Jt(z) / <Jo{z). 

Claim 5.4. \E(D(x, y, e)) - eE(Sf X! yf)\ < ||/||oo(6eP(E) + e 2 ). 


Proof. This is claim is justified by playing around with definitions, but let us do it in 
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detail nonetheless. Let E\ c E be the event that for some z G [x,y\ and some t G [0, e], 
cTt(z) Y ctq(z), and in addition there was no Poisson arrival at x during the time interval 
[0, t]. Then 

E(D • l{£7i}) < 2||/|| 00 eP(£? 1 ) < 2||/|| 00 eP(£). ( 22 ) 

The first inequality is due to the fact that after E\ happens, it is still necessary to have 
a Poisson arrival in [i, e], which has probability less than e, independently of E\. Even if 
these events both occur, the maximum effect on the value of / is 211 /Hoq. 

Another case which is easily dispensed with is the event that there are two or more 
arrivals at x (denote the number of arrivals by k and this event by E 2 , so E 2 = {k > 2 }). 
Then IP(£ 2 ) < eP (E) because after the first arrival (which changes a x hence is included 
in E ) we still need another arrival, independently. Hence E (D ■ \{E- 2 }) < 2||/|| 00 eP(£). 
For the remainder (denote it by £3 = E fl (E\ U E 2 Y ), let 

G = (J? x , y f)(a 0 )l{k = l}. 


Then 

E(D • 1 {£ 3 }) = E(G • 1 {£ 3 }) 

because under £3 these are exactly the same variables. Hence 

|E(£• 1{£ 3 }) — E(G)| < 2II/IUP ({k = 1}\£ 3 ) = 2||/|| 00 P({fe = l}n£i) < 2||/|| 00 eP(£) 

where the equality follows because {k = 1} C E\E2, and the last inequality is as in (22). 
Finally the definition of Jz f x ,y gives E(G) = ee” e E E£ x<y f so |E(G) — eEJ£ Xjy f\ < e 2 (recall 
that we assumed e < ^). Putting everything together the claim is proved. □ 

Differentiability of E D(x, y, t ) is now immediate; we write 

y— 1 

E = { > 0} 

z=x 

and get from Lemma 5.2 and Markov’s inequality that P(£) < eC(x,y). Claim 5.4 now 
gives that ^E D(x,y,t) = Ej V x<y f. 

To be able to sum the derivatives over x we use the assumption that E(^(cro) 1+?7 ) < 00 
from the local Condition A and Markov’s inequality. We conclude that for any 5 > 0 and 
j/GZ there exists N = N(y ) such that 

E (l y (a t )l{ly(a t ) > N}) <5 Vi G [0,i 0 ]. (23) 
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Assume also N > i. Applying claim 5.3 to all x < y — N and summing gives 

^2 ^D(x,y,e) < 2||/||oo f ly{a t )l{ly(a t ) > N} dt < 2||/||oo • eS 
x<y—N ^ 

SO 

e—>o e E 2/? ^ ^ ^^x,yf H _ 2||/||oo5 

X <V y,x<y—N(y) 

lim-V ] D (x,y,e)> V E^/-2||/|| 0O 5. 

f—£ — 

Taking A —>• 0 gives the lemma. □ 

Remark 5.5. The proof of (4) just given can be strengthened slightly when process is 
stationary. In this case it is enough to assume E(7 x (cr(0))) < C(x) for all x i.e. it is not 
necessary to have 1 + 7/ moments, 1 is enough. The only difference in the proof is the 
justification of (23). 

Let us now give the proof of the conditional version of the last lemma. 

Proof of Lemma 2.11, (3). Let x < y with y in the support of /, and let e € (0, |). We 
wish to define the effect of Toorn jumps from x to y on f(a e ). Therefore let t\,..., tk be 
the times of the Poisson arrivals at x in the time interval [0, e] and define, for x < y, 

k 

D(x,y,e ) = ^(/(<r tj ) - f(a t -))t{a s (x) = • • • = a s {y - 1) + er s (y)}. 
i =1 

Since, for all z, a Toom process changes at z only finitely many times in any time interval 
and / is local, we get f(a e ) — f(a o) = Yh x<y D{x, y, e). Let f(a) be the event that 
a(x) = a(x + 1) = ••• = a(y — 1). We define two “bad” events, B 1 and B 2 as follows. 

1. For z £ [x,y] let B 1 (x, z,y,t) be the event and that for some s < t we have 
a s (z) f &o(z); and that for some u € (s,t) there was a Poisson arrival at x and 
£(<7 u -) occurred. 

2. Let B 2 (x,y,t) be the event that £(o y) occurred, that for some s < t. we have 
cr s (x) o’o(x), and that there is a Poisson arrival at x in the time interval (, s,t ). 
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Let 


B l [t) = |J B 1 (x, z, y, t) 

x<z<y: 
yeSupp(f) 

B 2 [t) = U B 2 (x, y, t) 

x<y: 

yeSupp(f) 

Claim 5.6. Almost surely, 

sup \E(f(a t ) - /(cj 0 ) | cr 0 ) - t(5?f)(ao)\ 
te[o,e] 

< 2||/|| 00 (P(S(e) | cr 0 ) + e 2 )^l + ly( a o))- 

2 /eSupp(/) 

Proof. Fix one y in the support of / and one x < y. Define the variable G(x, y, t ) = 
{J£ x . y f){cr o) • 1 {k > 1} where k is the number of Poisson arrivals at x in the time interval 
[0, t\ ( G depends on t. only via k). We claim that if B(e ) doesn’t happen then D(x , y, t ) = 
G(x, y, t ) for all t € [0, e]. 

Indeed, if k = 0 then they are both 0, and if k > 1 then G = (J'? X: yf)(cro) while 
D = ^T(2zf Xiy /)(<r,-), w h ere > as usual, ti are the Poisson arrivals at x in the time interval 
[0, t]. If B 1 doesn’t occur, only the first term in this sum may be non-zero (this first 
arrival will change the value of er(x), so any further arrival at x, if it contributes to D it 
must also trigger B 1 (x, x, y, e)). Further, the fact that B 1 doesn’t occur implies that 

0 =S> a t^\ [xy] = a v\[x,y\ 

while the fact that B 2 doesn’t occur implies that 

(J? x , y f)(a 0 )^0 => a t-\ [xy] = a v\[x,y\ ■ 

Thus if neither occurred, = {J£ x ^ y f){op) and D = G. 

Summing over x and y gives, assuming B = B(e) doesn’t occur, that 

f{°t) - f{o o) = G{t) := ^2 G ( x ’ 2b 0 

x,y 

(note that the sum defining G is in fact finite). Taking conditional expectation gives 

mn°t) - /(<7 0 )) • 1 {B c } I do) = E(G • 1 {B c } I <t 0 ) (24) 


B(t) = B 1 {t)AB 2 {t). 
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One remainder can be estimated simply by 


E(\f(a t ) - /(<t 0 )| • 1 {B} | <T 0 ) < 2\\f\\ 00 F(B \ a 0 ) (25) 

while for the other we ignore the condition {k > 1} in the definition of G and get 

E(|G| • 1 {B} | <X 0 ) < e( Y \& x ,yf((To)\ ■ 1 {B} | <7 0 ) 

X,y 

= Y \&x,yf{oo)\ • HB | <tq) < 2H/IU Y lyfroWB I <70) (26) 
x,y y 

Summing (24), (25) and (26) gives 

|E(/fo) - f(<T o) I <7 0 ) - E(G Ml < 2||/|| 00 P(5 I cr 0 )(l + Y l v( a o)) 

y 

Finally a Poisson process calculation gives E(G(x, y, t) | cro) = (1 — e~ t )^f x ^ y f(ao) so 
|E(G|cr 0 ) — t(3f f)(ao )I < 2 ||/||ooe 2 lOy ^(<7o)- Putting everything together the claim 
is proved. □ 

Thus the lemma will be proved once we estimate P (B). We start with B 1 . 

Claim 5.7. With probability 1, 


lim 2 n P(L> 1 (2 _n ) | cj 0 ) = 0. 

n—>oo 

Proof. Fix x < z < y. We note two estimates for B l (x,z,y,e). First, ignoring the 
requirement at z gives 

PMOo z, Vi e)) < e sup P (l y > y - x) < C(y) ■ e(y - x)^~ v 

t£ [0,e] 

where the second inequality is due to our moment assumption on a and Markov’s in¬ 
equality. Second, ignoring the requirement that the arrival at x actually changes y, and 
using only the fact that an arrival happened we get from Lemma 5.2 

nB\x,z,y,e))<C{z)-e 2 . (27) 


Summing over all x gives 

p(U ^(abz-y-e)) <e^min (C(y)(y - x) -1-11 ,C(z)e) < C{y, z)e 1+v/(1+v) 

X X 

where C(y,z) is some constant which depends on y and z but not on e. Let us denote 
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B 1 (z,y,e ) = U x B 1 (x, z,y, e). Markov’s inequality now gives for the conditioned events 
that 

¥{¥(B 1 (z,y,e)\a 0 ) > < C(y, z)e Tl ^ 1+ri \ 

This implies, by the Borel-Cantelli Lemma, that with probability 1, 

lim 2 n P(B\z iyi 2- n )\a 0 ) = 0 (28) 


for all z and y. 

There is another important consequence of the conditioning over oq- Fix some value 
of do and let z* = z*(y) = y — l y (a o) — 1. For any x < z* we have that 

y 

1 J B 1 (x,z,y,e) C B 1 (x, z*,y, e) 

Z=X 

because the requirement that at the time u (u from the definition of B 1 ) we have the 
event that £(cr u -) occurs is fulfilled only if cr u -(z*) / ao(z*). Thus 


lim 2 n p( M M B\x,z,y,2- n ) 

n—>-oo \ v - / 

x<z* z 


do < lim 2 "P( j B 1 (^, 2/ ,2-)|(7o) = 0. 
/ n—> oo 


The remaining x (i.e. x > z*) can be estimated directly by summing (28) over z from 
z* + 1 to y. We get 


I/O— n\ 


B\2 


^0 


<2 n Y, E nBHz,y,2- n )\a o )^0 

yGSupp f z=z* (y) 


proving the claim. 


□ 


The estimate of B 2 (e) is much simpler and we will not dignify it with a claim. By 
Lemma 5.2, 

nB 2 \a 0 )<e 2 ]T C(x) 

2/eSupp / x=y-l y (<T 0 ) 


and hence 


lim 2 n P(5 2 (2" n ) | (j 0 = 0 

n—> oo 


(29) 


almost surely. 
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The lemma is now proved. Applying claim 5.6 for t £ [2 n ,2 "] we get 


- f (cr 0 ) | cr 0 ) - - 27 (cr 0 ) 

< 2||/|| 00 (2 n+1 P(B(2- n ) | <t 0 ) + 2 _n+1 ) (l + £ Z y (u 0 )) . 

y 

By claim 5.7 and (29), 2 n ¥(B(2~ n ) \ cfq) —> 0 almost surely, and the other terms on the 
right hand side are constant for any fixed <jq. □ 
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